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Heavy-light meson in anisotropic lattice QCD* 
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We examine whether the 0(a) improved quark action on anisotropic lattices can be used as a framework for 
the heavy quark, which enables precision computation of matrix elements of heavy-light mesons. To this end, it 
is crucial to verify that a mass independent and nonperturbative tuning of the parameters is possible. As a first 
step, we observe the dispersion relation of heavy-light mesons on a quenched lattice using the action which is 
nonperturbatively tuned only for the leading terms. On a lattice with the spatial cutoff ^ — 1.6 GeV and the 
anisotropy ^ = 4, the relativity relation holds within 2% accuracy in the quark mass region aaVriQ < 1.2 with the 
bare anisotropy parameter tuned for the massless quark. We also apply the action to a calculation of heavy-light 
decay constants in the charm quark mass region. 
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1. Introduction 

Recent experimental progress at B factories 
suggests that precise computation of hadronic 
matrix elements in lattice QCD is a key to the 
search for signals of new physics in flavor physics. 

However, the two major formalisms for the 
heavy quark, i.e. HQET approaches or the rela- 
tivistic approaches to describe the heavy quarks 
on the lattice still suffer from perturbative and/or 
discretization errors, which are typically ~ 10%. 
We therefore need yet another formalism of the 
heavy quark for the precise and systematic com- 
putation of matrix elements. 

Ideally, in such a formalism, one should be able 
to (i) take the continuum limit, (ii) compute the 
parameters in the action and the operators non- 
perturbatively, (iii) and compute matrix elements 
with a modest computational cost. 

As a candidate of framework which fulfills these 
conditions, we investigate the anisotropic lattice 
with the temporal lattice spacing Ur finer than 
the spatial one Oo- Anisotropic lattice 

framework evidently satisfies above conditions (i) 
and (iii). Our expectation is that on anisotropic 
lattices the mass dependence of the parameters 
becomes so mild that one can adopt coefficients 
determined nonperturbatively at massless limit. 
Whether these promises will be practically satis- 
fied should be examined numerically. 



As a first step toward this goal, we study the 
action which is nonperturbatively tuned only for 
the leading term, while 0(a) improvement terms 
are tuned at the tadpole tree level. We measure 
the heavy-light meson dispersion relation with the 
anisotropy parameter tuned at the massless limit. 
We can then study the breaking of relativity con- 
dition as a function of the heavy quark mass, with 
which we can test if the small mass dependence 
of the parameter really holds. 

2. Anisotropic lattice quark action 

Our heavy quark formalism basically follows 
the Fermilab approach M but is formulated on 
the anisotropic lattice 

The quark action is represented as 



Sf ^tlj{x)K{x,y)ilj{y) 
K{x,y) = Sx^y - K 



(1) 



il-j4)U4{ x)5^^^ y 

+ (l+74)C/i(2:-4)<5^_4,, 

l^<rT,i [ ir-l^)U^{x)S^^■. y 

+ {r + 7^)uUx-i)S.^^l^y 

KaCE}2i'^iiFii {x)6x,y 

rK„CBj2i>j(^i]Fij{x)dx^y, (2) 
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where Ko- and k^- are the spatial and temporal 
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hopping parameters, r is the Wilson parameter 
and Ce and cb are the clover coefficients. 

In principle for a given Ko-, the four parame- 
ters Ka/Kr, T, Ce and Cb should be tuned so that 
Lorentz symmetry holds up to discretization er- 
rors of O(a^). 

In this work, we set the spatial Wilson param- 
eter as r — 1/^ and the clover coefficients as 
the tadpole-improved tree-level values as ce — 
l/uaU^, Cb = l/i^CT ^'^d perform a nonperturba- 
tive calibration only for 7^. The tadpole improve- 
ment Q is achieved by rescaling the link variables 
as Ui{x) — > Ui{x)/ua and U4^{x) U4^{x)/ut, 
with the mean-field values of the spatial and tem- 
poral link variables, Mo- and u-n respectively. This 
is equivalent to redefining the hopping parame- 
ters with the tadpole-improved ones (with tilde) 
through /ua and Ht = kr/ur- We define 

the anisotropy parameter as 7f = Kt/k,(t- 

For convenience, we also introduce k as 

- = ^ -2(7F + 3r-4) =2(mo7F+4), (3) 

where toq is the bare quark mass in temporal 
lattice units. 7^^ is the bare anisotropy param- 
eter, which should be tuned so that the physical 
fermionic anisotropy S.F retains the same value 
as that of the gauge field, <^g- For the quenched 
case, as for the present work, one can first cali- 
brate the anisotropy parameter of gauge field sep- 
arately from the quark field, and then tune 7f so 
that ~ S.G ^ ^ holds. 

The calibration of bare anisotropy 7f was per- 
formed in Ref. |^ using the dispersion relation 
of mesons for quark masses less than the charm 
quark mass region, on quenched lattices with 
^ = 4 for three different lattice spacings. The 
quark mass dependence is indeed mild and well 
represented by a linear form in m^. 

3. Relativity condition 

We next study whether the relativistic disper- 
sion relation holds for the heavy-light mesons 
with 7i? tuned for the massless quark In gen- 
eral, the lattice dispersion relation for arbitrary 
values of is described as 
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Figure 1. The heavy quark mass dependence 
of the anisotropy for heavy-light and heavy- 
heavy mesons at /3 = 5.95. 



where the energy E and the rest mass m are in 
temporal lattice units while the momentum p is in 
spatial units. The parameter in this equation 
defines the anisotropy of the quark fields, whose 
difference with the gauge field anisotropy ^ probes 
the breaking of relativity. 

Numerical simulations are performed on a 
quenched lattice of the size 16'^ x 128 generated 
with the Wilson plaquette action at /? = 6.10 
for ^ = 4, which is realized by taking the bare 
anisotropy of gauge field 70 using the numerical 
formula determined by Klassen in one percent ac- 
curacy . The lattice cutoff defined through the 
hadronic radius ro is = 1.623(9) GeV. Other 
details of parameters such as the mean-field val- 
ues are found in Refs. ||^,^. 

The result for the pseudoscalar meson is dis- 
played in Fig. |l|. Similar result is obtained for 
the vector channel. The measured value £,p for 
the heavy-light meson coincides with ^ in the 
quark mass region niq < 0.3 (in temporal lat- 
tice units) within 2% accuracy. (This accuracy is 
just a guide for argument of uncertainty in this 
paper.) Therefore, for the range of quark mass 
TTtq < 0.3 which covers the charm quark mass 
with present lattice cutoff, our quark action with 
the bare anisotropy tuned for massless quark can 
successfully describe the heavy-light meson with 
controlled systematic uncertainty. We also show 
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the of the heavy quarkonium. In this case, 
rapidly increases because of large momentum 
transfer inside heavy quarkonium, p ~ amq, in 
contrast to the case of heavy-light systems, where 

P — ^QCD- 

4. Heavy-light decay constant 

The result in the previous section tells us 
that the framework can be applied to charmed 
heavy-light meson systems. As a first appli- 
cation, we compute the heavy-light decay con- 
stants. Numerical simulations are performed on 
two quenched lattices at /3 = 5.95 and 6.10 with 
the renormalized anisotropy ^ = 4. The former 
is the same as in previous section, and the lat- 
ter has spatial cutoff a^^ = 2.030(13) GeV. For 
the light quark, we use three values of the hop- 
ping parameter k which cover the range of quark 
masses from m<; to 1.5 m^. The chiral limit is 
taken by a linear extrapolation. The matching 
to the continuum theory is at the tadpole im- 
proved tree level. The result at /3 = 6.0 is dis- 
played in Fig. Ij. In this analysis, the physical 
quark masses are defined through the lattice cut- 
off defined with the K* meson mass. We find 
that the behavior of fy/m is consistent with the 
results on isotropic lattices. Since the currents 
are matched only at the tadpole tree level and 
the 0(a) improvement terms are tuned also only 
at the tadpole tree level, the present calculation 
contains rather large renormalization uncertainty 
as well as the cutoff dependence. In order to sup- 
press these errors we compute the ratio of decay 
constants. Our preliminary results are 

fo/f. = 1.566(43) at /3 = 5.95, 

1.515(43) at /3 = 6.10, 

foJfD = 1.140(14) at /3 = 5.95, 

1.142(14) at /3 = 6.10, (5) 

which are consistent with previous works on 
isotropic lattices [^. A remarkable point is that 
the lattice spacing dependence of the ratio is 
small. In fact, individual decay constants could 
suffer from both the mass independent and mass 
dependent errors, but the former can largely can- 
cel in the ratio. Therefore, the stability of the 
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Figure 2. Heavy- light meson decay constant in 
the pseudoscalar channel at /3 = 6.10. 

ratio may be an indication that the mass depen- 
dent errors are suppressed in our formalism on 
anisotropic lattices as expected. 

In conclusion, the results of numerical simu- 
lations in quenched lattices are encouraging for 
further development in this direction. Further 
improvements, such as nonperturbative tuning of 
the clover coefficients, are necessary for achieving 
the desired accuracy. 

The simulations were done on a NEC SX-5 at 
RCNP, Osaka University and a Hitachi SR8000 
at KEK. 
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